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Horava and Melby- Thompson recently proposed a new version of the Horava-Lifshitz theory of 
gravity, in which the spin-0 graviton is eliminated by introducing a Newtonian pre-potential ip and 
a local (7(1) gauge field A. In this paper, we first derive the corresponding Hamiltonian, super- 
momentum constraints, the dynamical equations, and the equations for <p and A, in the presence of 
matter fields. Then, we apply the theory to cosmology, and obtain the modified Friedmann equation 
and the conservation law of energy, in addition to the equations for if and A. When the spatial 
curvature is different from zero, terms behaving like dark radiation and stiff-fluid exist, from which, 
among other possibilities, bouncing universe can be constructed. We also study linear perturbations 
of the FRW universe with any given spatial curvature k, and derive the most general formulas for 
scalar perturbations. The vector and tensor perturbations are the same as those recently given by 
one of the present authors [A. Wang, Phys. Rev. D82, 124063 (2010)] in the setup of Sotiriou, Visser 
and Weinfurtner. Applying these formulas to the Minkowski background, we have shown explicitly 
that the scalar and vector perturbations of the metric indeed vanish, and the only remaining modes 
are the massless spin-2 gravitons. 

PACS numbers: 04.60.-m; 98.80.Cq; 98.80.-k; 98.80.Bp 



I. INTRODUCTION 

Recently, Horava proposed a quantum gravity theory 
1], motivated by the Lifshitz theory in solid state physics 
2| , for which the theory is often referred to as the Horava- 
Lifshitz (HL) theory. It is non-relativistic and power- 
counting ultraviolet (UV)-renormalizable, and was ex- 
pected to recover general relativity (GR) in the infrared 
(IR) limit. HL theory has attracted a great deal of at- 
tention due to its several remarkable features, such as 
the divergence of its effective speed of light in the UV, 
which could potentially resolve the horizon problem with- 
out invoking inflation [3|. Scale- invariant super-horizon 
curvature perturbations can also be produced without in- 
flation pl-tlllj , and dark matter and dark energy can have 
their geometric origins [l2l [l3j. Furthermore, bouncing 
universe can be easily constructed due to the high-order 
derivative terms of the spacetime curvature |14l - [la | . For 
detail, we refer readers to [l7| and references therein. 

The HL theory is based on the perspective that Lorentz 
symmetry should appear as an emergent symmetry at 
long distances, but can be fundamentally absent at high 
energies [llj]. With this in mind, Horava considered sys- 
tems whose scaling at short distances exhibits a strong 
anisotropy between space and time, 



£x, t -> £ z t. 



(1.1) 



In (d + l)-dimensional spacetimes, in order for the the- 
ory to be power-counting renormalizable, it requires z > 
d. At low energies, the theory is expected to flow to 
z = 1, whereby the Lorentz invariance is "accidentally 
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restored." Such an anisotropy between time and space 
can be easily realized when one writes the metric in the 
Arnowitt-Deser-Misner (ADM) form [l]|, 



ds 2 



-N 2 c 2 dt 2 + gij (dx l + N l dt) (dx j + N J dt) , 

(t, j = 1,2,3). (1.2) 



Under the rescaling (jl.ll) with z = d — 3, a condition 
we shall assume in the rest of this paper, the dynamical 
variables 7Y, N l and gij scale, respectively, as, 



N —> N, N l ^£- 2 N\ 9ij > 9ij . 



(1.3) 



The gauge symmetry of the system are the foliation- 
preserving diffeomorphisms Diff(Af, J 7 ), 

t = t-f(t), £ 4 =.T l -0(i,x), (1.4) 
for which the dynamical variables change as 
Sgij = ViCj + V,Ci + f9ij, 



SNi = A^V 4 C fc + ( k V k N t + g lk ( k + NJ + Nj, 



SN = ( k V k N + Nf + Nf, 



(1.5) 



where / = df/dt, V; denotes the covariant derivative 
with respect to the 3-metric g^, and Ni = gi k N k , etc. 
From these expressions one can see that the lapse func- 
tion N and the shift vector Ni play the role of gauge fields 
of the Diff(M, F) symmetry. Therefore, it is natural to 
assume that and Ni inherit the same dependence on 
spacetime as the corresponding generators, in addition 
to the fact that the dynamical variables g^ should in 
general depend on both time and space, that is, 

N = N(t), Ni = Ni(t,x), <i,, //,,;/../•). (1.6) 

which is clearly preserved by the Diff(AT, J 7 ), and often 
referred to as the projectability condition. 
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Due to the restricted diffeomorphisms (jl.4l) . one more 
degree of freedom appears in the gravitational sector - 
a spin-0 graviton. This is potentially dangerous, and 
needs to decouple in the IR regime, in order to be con- 
sistent with observations. Unfortunately, it was shown 
that this might not be the case. In particular, the spin- 
mode is not stable in the original version of the HL 
theory [l| as well as in the Sotiriou, Visser and Weinfurt- 
ner (SVW) generalization 0, H[. Note that in both of 
these two versions, it was all assumed the projectability 
condition. In addition, these instabilities are all found 
in the Minkowski background. Recently, it was found 
that the de Sitter spacetime is stable in the SVW setup 
[l7j . So, one can take the latter as its legitimate back- 
ground, similar to what happened in the massive gravity 
2C|. However, the strong coupling problem still exists 
21 , [22 1 , although it might be circumvented by the Vain- 
shtein mechanism [23j , as recently showed in the spherical 
[HJ and cosmological [13] cases. 

On the other hand, giving up the projectability condi- 
tion, that is, assuming that the lapse function N depends 
on both time and spatial coordinates, Bias, Pujolas and 
Sibiryakov (BPS) [25[ found that inclusion of terms made 
of a,-, 



take exactly its relativistic value A = 1. Note that in GR 
the spacetime diffcomorphism symmetry, Diff(Af), 



Oi = dAn(N), 



(1.7) 



can cure the instability of the Minkowski spacetime. By 
properly choosing the coupling constants, the strong cou- 
pling problem (26l - [29j can be also addressed [30j . How- 
ever, a price to pay is the enormous number of indepen- 
dent coupling constants: only the sixth-order derivative 
terms in the potential are more than 60 [29]. It should 
be also noted that giving up the projectability condition 
often causes the theory to suffer the inconsistence prob- 
lem [3lj | . Kluson recently showed that the Hamiltonian 
formalism of the BPS model is very rich, and the corre- 
sponding algebra of constraints is well-defined [32| . 

To cure the instability problem, another very attractive 
way is to eliminate the spin-0 graviton from the theory, so 
that the resulting one has as many generators per space- 
time point as GR does. This is done recently by Horava 
and Melby- Thompson (HMT) [H| (with the assumption 
of the projectability condition (|1.6I) ) by extending the 
foliation-preserving-diffeomorphisms, Diff(Af, 7), to in- 
clude a local £7(1) symmetry, 



17(1) ix Diff(M, 7"). 



(1.8) 



Effectively, the spatial diffeomorphism symmetries of GR 
are kept intact, but its time reparametrization symmetry 
is linearized and the corresponding algebra is contracted 
to a local gauge symmetry [34j . The restoration of gen- 
eral covariance, characterized by Eq. (jl.8p . nicely main- 
tains the special status of time, so that the anisotropic 
scaling (jl.ip with z > 1 can still be realized. 

A remarkable by-production of this "non-relativistic 
general covariant" setup is that it forces the coupling 
constant A, introduced originally to characterize the de- 
viation of the kinetic part of the action from GR [l[ , to 



jgM=xM x)j (^ = 0,1,2,3) 



(1.9) 



also forces A = 1 and protects this value from quantum 
corrections. 

At short distances, the theory exhibits a high 
anisotropy between time and space. As a result, the UV 
behavior of the theory is dramatically imporved. At long 
distances, the theory is driven to an IR regime, where it 
shares many features with GR. In particular, under the 
influences of the relevant terms, the scaling is naturally 
isotropic with the relativistic value z = 1. Moreover, 
since the extended symmetry forces A = 1, in the IR limit 
the action will be dominated exactly by the Einstcin- 
Hilbert terms in the ADM decomposition [19j . 

In this paper, we investigate this new version of the HL 
theory. Specifically, in Sec. II we first give a brief review 
of it, and then derive the corresponding Hamiltonian, 
super-momentum constraints, the dynamical equations, 
and the held equations for the Newtonian pre-potential 
<p and the U(l) gauge field A, in the presence of mat- 
ter fields. The potential used in this paper is the one 
constructed by SVW Q > which represents the most gen- 
eral form, subjected to the assumptions that it respects 
the parity and its highest order of the spatial derivatives 
is six, the minimal requirement to have the theory be 
power-counting renormalizable [l|. In Sec. Ill, we ap- 
ply the theory to cosmology, and obtain the modified 
Friedmann equation and conservation law of energy, in 
addition to the equations for (p and A. In Sec. IV we 
study the linear perturbations of the FRW universe with 
any given spatial curvature k, and present the general 
formulas for scalar perturbations. The vector and tensor 
perturbations are the same as those given by one of the 
current authors in [ll| in the SVW setup 0] , because the 
gauge field and the Newtonian pre-potential have no con- 
tributions to these parts. Applying these formulas to the 
Minkowski background in Sec. V, we show explicitly that 
the scalar and vector perturbations of the metric vanish 
identically. The only non-vanishing dynamical variables 
are the traceless and divergence-free tensor H%j, which 
describes the massless spin-2 graviton, a situation that is 
precisely the same as in GR. These results are consistent 
with the ones obtained earlier in (33|. In Sec. VI, we 
present our main conclusions. 



II. NON-RELATIVISITC GENERAL 
COVARIANT HL THEORY 

In order to limit the spin-0 graviton, HMT introduced 
two new fields, the f7(l) gauge field A and the Newtonian 
pre-potential tp, where in general both of them depend 
on space and time, 



A = A(t,x k ) 



(2.1) 
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Note that the notations used in th is p aper are slightly 
different from those adopted in [33| 1 . Under the 
Diff(Af, J 7 ), these fields transfer as, 



5 A = (%A + fA + fA, 
dip = f<p + Cdi<P, 



(2.2) 



while under the local U(l), they, together with g. t j, trans- 
fer as 

S a A = a — AT'ViQi, 5 a ip = —a, 
5 a Ni = NVia, 8 a9l] = = S a N, (2.3) 

where a is the generator of the local U (1) gauge symme- 
try. For the detail, we refer readers to [33j. 
The total action is given by, 



S 



dtd A x 



Ck — C\ 



C v + Ca 



(2.4) 



where g = det gij , and 



C K = K t jK %3 - K 2 



C v = ygv^Kij + ViVjtp 



(2.5) 



Here the coupling constant A 9 , acting like a 3- 
dimensional cosmological constant, has the dimension of 
(length) ~ 2 . The Ricci and Riemann terms all refer to the 
three-metric g^ , is the extrinsic curvature, and Qij is 
the 3-dimensional "generalized" Einstein tensor, defined, 
respectively, by 



of the coupling constants. With this condition, Cy takes 
the simple form, 

C v =w 2 C ij C ii , (2.8) 

where w is a coupling constant, and Cy denotes the Cot- 
ton tensor, defined by 



(2.9) 



In 0, by assuming that the highest order derivatives 
are six and that the theory respects the parity, SVW 
constructed the most general form of Cy , given by 

Cv = C 2 9a+9iR+^(92R 2 +g 3 R lJ R 13 ) 

+^ (94R 3 + g 5 R /.'... /.'" + nuirj!! !! 1 ; 
+i [gyRy 2 R + gs (V 4 i? jfc ) (V 4 i?^)] , (2.10) 

where the coupling constants g s (s = 0, 1, 2, . . . 8) are all 
dimensionless. The relativistic limit in the IR requires 
gi = —1 and £ 2 = l/(167rG) [7[. In this paper, we 
shall be concerned only with this potential, and our for- 
mulas to be obtained below can be easily generalized to 
other forms of the potential, including the one given by 
Eq.flUS]), and the f(R) term studied in pj]. 

Variation with respect to the lapse function N(t) yields 
the Hamiltonian constraint, 

J d 3 x^/g~ (C K + Cy - tpG ij ViVj<p) =8ttG J d 3 x^fg~j\ 

(2.11) 

where 



K, 



1 

2N 



i-9ii 
1 



V 4 TV, 



:9ijR + A g9ij- 



(2.6) 



Cm is the matter Lagrangian density, which in general 
is a function of all the dynamical variables, U(l) gauge 
field, and the Newtonian prepotential, 



C M =C m (N, Ni, g^, ip, A; x), 



(2.7) 



where x denotes collectively the matter fields. Cy is an 
arbitrary Diff(I])-invariant local scalar functional built 
out of the spatial metric, its Riemann tensor and spatial 
covariant derivatives, without the use of time derivatives. 
In the original approach of Horava [lj, the detailed bal- 
ance condition was imposed, in order to limit the number 



jt = 2 HNC M ) 



SN 



(2.12) 



Variation with respect to the shift N l yields the super- 
momentum constraint, 

Vj (n ij - tpg ij ^j = 8irGJ\ (2.13) 

where the super-momentum -k 13 and matter current J 1 
are defined as 



7r ij = S ( NCkS) = + Kg iS , 



J' 



Sgij 
5C M 



-N 



SN,, ' 



(2.14) 



Similarly, variations of the action with respect to ip and 
A yield, 



In particular, we have ip 



Ki 



Q HMT , Qij = S^j MT , where quantities with super indice 
"HMT" are the ones used in [H. 



G i] [h,, ■ V,V,,;J = 8ttGJ v , 
R=2A q + 8irGJ A , 



(2.15) 
(2.16) 
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where 



T SC M t _J(NC m ) 



RfRlRZ-SVkVt {R 2 f 



Sip 



6A 



(2.17) {F 7 ) l0 = 2V 4 V, (V 2 i?) - 2 (V 2 i?) R^ 



On the other hand, variation with respect to gij leads to 
the dynamical equations, 



— yg (*« - ^«)] t = 2 ( K y + 2KK* 



H Vfe 

JV 

1 



iV fc 7r ij ' - 2n k ^N j ^ 



+ -{C K + C v + C A )g t3 
+ F« + F« + F£ + 87rGr y , (2.18) 

where (if 2 ) y = A**/, = (/,, + f Jt ) /2, and 



s=0 



V5 <%j 

3 

— \ ^ 

V ~ L^i (¥>,")' 
n=l 

= [AR ij - -g ij V 2 ^A 



(2.19) 



with n s = (2, 0, -2, -2, -4, -4, -4, -4, -4). The stress 
3-tensor t 4 - 7 is defined as 



r « = 2 <5(VI7£m) ; 
V5 <faij 



(2.20) 



and the geometric 3-tensors (-F s )jj and FV > are defined 
as, 



(*b)« 



1 



(F 2 ) tf - 2 (i? 4 , - ViV,-) R - -9a (R - 4V 2 ) R, 
(Fshj = V 2 i? y - - (V 4 V, - 2,R l3 ) R 4 (R 2 ).. 

+ (3i? fe /i?' £ ' + V 2 i? - 2i? 2 ) , 

(fkjy = 3(R ij -V i V j )R 2 -±g ij (R-6V 2 )R 2 , 
- ( J R y +V l V J )( J R feJ i? fe ')+2i?(i? 2 ) jj 

+V 2 (iifly) - V fe [V, {RR jk ) + Vj (ilflifc)] 

-i 5ij [(i?-2V 2 ) (JfcJl«) 

-2V fc V, , 
= 3 (*% + | [v 2 (i? 2 V . 



(V,i?) (V,i?) - [(VR) 2 + 4V 4 i? 



(^8)y 



T?i3 



= V 4 i? 2J - V fc (V,V 2 i?* + Vj-V 2 i2?) 
-(ViiZf) (vX) - 2 (v fe i?0 (V fc %) 
-\9ij [(V fc i? im ) 2 - 2 (VfeViV 2 ^)] , (2.21) 

= \p { (W + V 2 p) i# - 2 (2^ + V J V fe ^) i? lfe 
-(2A g -R) (2^'+V l VV)}, 

^v fe {^(^ + vv) 

+ ^(^ + w)-^(^ + w)}, 

= 1 {2V fe V<7£ )fe - V 2 .^' - (VfcVj/*') g 4 -' } , 

(2.22) 



^(V,3) 



where 



/« = ^ I (2^ + v*v^) - i(2A- + vvy 

(2.23) 

The matter quantities (J*, J 2 , J v , Ja, t u ) satisfy 
the conservation laws, 



1 



A 



2^fe , r- Tk , 



N^g 



= 0, (2.24) 



-^V fc J fc + J v V^-^V^-0. (2.25) 



III. COSMOLOGICAL MODELS 



The homogeneous and isotropic universe is described 



by, 



where 



N = l, N l= 0, g ij =a 2 (t) lij , (3.1) 



7y (i+^r 



(3-2) 



-v fe (V, (i? 2 ) jfe 



with r 2 e i 2 + 1/ 2 + z 2 and fc = 0,±1. Using the U(l) 
gauge, on the other hand, we can set 



<p = 0, 



(3.3) 
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without loss of generality. Then, we find that 



Ki. 



-a Hj i: 



Ri 



(3.4) 



where H = a/a. Thus, we obtain 
C K = -QH 2 , C v = 0, 

£k = 2A+ ^ + 1?^! + ^, (3. 5) 
a z or a° 

where A = ( 2 g /2, and 

a 3g 2 +g 3 a 9.94 + 3.9s + .96 ,„ «\ 

P\ = — ^ , P2 = ^4 • (3.6) 

The matter components are 

J* = -2p, J 4 = 0, T ij =pg ij , (3.7) 

where p and p are the total density and pressure of the 
matter fields. Then the Hamiltonian constraint (|2.11[) 
reduces to the super-Hamiltonian constraint, 

C K (t)+C v {t) = 87rGJ*(t), 

which leads to the modified Friedmann equation, 

2 gik _ 8ttG A 2/3ifc 2 4/3 2 fc 3 



3^3 



From Eqs. ([2TT4]) and we also find that 



(3.8) 



7r« = -2Hg ij . (3.9) 

Hence, the dynamical equation (|2.18l) reduces to Q 

4ttG , „ s 1 . A / . fc 

^2 



a 
n 



-( P + 3p) + -A--(A 9 -- 



2/3ifc 2 8/3 2 fc 3 



(3.10) 



Similar to GR, the super-momentum constraint (I2.13[) 
is then satisfied identically, since J 1 = and, from 
Eq. p.6[) , VjTr 13 = 7r y ij = 0, where V, denotes the co- 
variant derivative with respect to 7^. Using Eqs. p.8[) 
and p.lOj) . it follows that in the FRW background the 
matter satisfies the conservation law, 



p + 3H (p + p) = AJ^ 



(3.11) 



Thus, due to the interaction between the gauge field and 
the fluid, its energy in general is not conserved. 

On the other hand, Eqs. (12.151) and (|2.16p yield, respec- 
tively, 



rr 1 • k \ 8lrG T 
//|A «-^J =- — J ^ 



3A- 



A„ = AuGJa- 



(3.12) 
(3.13) 



When matter is not present, we have J v = = J a- 
Then, Eqs. ([3-T2]l and (|3~T3]) implies that A s = = k, 
while Eq.([33J yields a(t) = e Ht , where H = ^AA 
which is the de Sitter spacetime. It is interesting to note 
that the de Sitter space can be also obtained from p = 
p = k = and J A = -A g /(4nG), J v = -3HA g /{8nG). 



It should be noted that the energy conservation p. lip 
can be also obtained from Eq. (|2.24j) . while the conserva- 
tion law of momentum, Eq. (|2.25l) , is satisfied identically. 
When /?i/?2 7^ 0, the corresponding terms act like a dark 
radiation and a stiff-fluid, respectively. Due to the pres- 
ence of these terms, one can easily construct bouncing 
universe in the early epoch of the universe EMI. 

In addition, in deriving Eq. p.8[) we followed the usual 
assumption that the whole FRW universe is homoge- 
neous and isotropic. In pjj], it was argued that such 
an assumption might be too strong. If one relaxes the 
assumption and requires that only the observed patch 
of our universe is homogeneous and isotropic, one can 
introduce the notion of "dark matter as an integration 
constant" of the Hamiltonian constraint (|2.11l) : p(t) in 
Eqs. (13.81) and (|3.11|i can be replaced by p(t)+£(t) in the 
observable patch, where £(t) = const /a 3 in the IR limit 
[l~2l . |H[ . Beyond the observable patch, £ is necessarily 
inhomogeneous. In order to analyze perturbations on an 
FRW background, one needs to restrict the perturbations 
to the observable patch, which then raises issues about 
matching across the boundary of the observable patch. 
In our approach, the background is a homogeneous FRW 
spacetime, so that £ = in the background. 



IV. COSMOLOGICAL SCALAR 
PERTURBATIONS 

In this section, we consider linear scalar perturbation 
of the FRW universe studied in the last section. We shall 
closely follow Q and use the notations adopted there 
without further explanations. However, in order to have 
the present paper as independent as possible, it is difficult 
to avoid repeating the same materials, although we shall 
try to limit it to its minimum. 

In the quasi-longitudinal gauge [8] , 

= = E, (4.1) 

the metric scalar perturbations are given by 

ds 2 = a 2 [-drj 2 + 2B\ i dx i dr] + (1 - 2tp) . 

(4.2) 

The gauge-invaraint quantities are now given by 

$ = HB + B', (4.3) 

where H. = a' /a and a prime denotes the ordinary deriva- 
tive with respect to 77. Using the U(l) gauge, we can 
further set 



S<P = 0, 



(4.4) 
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so that C v = 0. Following [8], we use quantities with 
over-bars as the ones calculated in the background. Then 
we find that 



Kij = 

Rij — 

Ck = 

C A = 

C v = 



-aH^ij + a [B {ij + (V/ + 2Hi/>) 7y] , 
6H 2 AH 



2A / 

A 9- 72 ) SA ~ — (vV + 3fc^ 



.91 



3jfe 
a 

4/3ifc 



-3k) if) 



3k) ip 



A8/3 2 k 2 



(4.5) 



where £y denotes the potential of the background given 
by Eq. (|3.5p . and B^ = V,, with Vj being the covariant 

derivative with respect to 7.^ and V 2 = 7 t5 V<Vj. 

To first-order the Hamiltonian constraint (|2.11[) takes 
the form, 



V7"<i 3 x 



V 2 + 3fc V " 



(2 - 3£)n 



(v 2 B + 3V>') 



6/3 2 & , 3.97 



C 4 a 4 
= 0, 



3k) ip 



(4.6) 



which is the same as that given in the SVW setup Q, 
where = —6J t /2. Eq. (|4.6p represents a generalization 
of the Poisson equation of GR [36] . 

To the first-order the supermomcntum constraint 
(|2.13[) takes the form, 



2ip' - 2kB = 8irGaq , 



(4.7) 



which is the generalization of the GR Oi constraint [36j, 



where SJ l = a 2 q' 1 . 



On the other hand, the linearized equations (|2.15l) and 
(|2~T6)) for the Newtonian pre-potential and the gauge field 
reduce, respectively, to, 



A,. 



a 

m 

a 2 



V 2 B + 3{i)' + 2Hip) 



VV + 3(2fc-a 2 A s )V> =8-kGcl5J ip , (4. 



\7 2 ip + 3kip = 2TrGa 2 5J A - 



(4.9) 



The linearly perturbed dynamical equations require 
the calculations of the perturbed (Fg)^ of Eq. (|2.2ip . 
which were given by Eq.(Al) in @. To avoid repeating, 
we shall not write them down here, and refer readers di- 
rectly to that paper. Then, we find that the trace part 



is given by 

i'' + 2Hip' - Ftp - ^l %3 5F l3 + - (v 2 B' + 2UV 2 B) 
+ 7^(v 2 -3A g a 2 + 6fc)v 

- g^(2V 2 + 3A s a 2 - 3kjSA = AnGa 2 5V, (4.10) 



where SF^ = J29sC ns ^( F s)ij, with £(F s )y given by 
Eq. (Al) in [|, and 



T = a 



-A 



g x k 2/3ifc 2 12/3 2 /c 3 



{5V + 2piP) j lj + U l{lj) 



(4.11) 



where the angled brackets on indices define the trace-free 
part: 



1 



l« 3 

The trace-free part of the dynamical equations is 



(4.12) 



B 



l(y> 



2HB 



1 



\m 



5F {ij) --(SA-A^) 
-8?rGa 2 n 



(ij) 



(4.13) 



Eqs. (|4.10p and (14. 13)) generalize the GR ij perturbed field 
equations [36]. 

The perturbed parts of the conservation laws Q2.24p 
and (|2~2"5|) give 



yfydPxl V + 3U (6T + 5(i)-3(p + p) tp' 



1 

2o^ 



(a 3 J a) 'SA + A^a 3 (S A - 3Atp)) ' 



= 0, 



(4.14) 



q' + 3Uq -aSV-^ (v 2 + 3fc) n + ^J A $A = 0, 

(4.15) 

where J a and J v are given by Eqs. (13.121) and (|3.13p . The 
energy conservation equation is an integrated general- 
ization of the GR energy equation, and the momentum 
equation generalizes the GR momentum equation 36] . 



V. STABILITY OF THE MINKOWSKI 
SPACETIME 

It can be shown that the Minkowski spacetime, 

a = l, A = if = k = 0, (5.1) 
is a solution of the HMT theory, provided that 

A g = A = J A = J v = p = p = 0. (5.2) 
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Then, the linearized Hamiltonian constraint (I4.6|) and the 
field equation (I4.8[) for Sip are satisfied identically, while 
Eqs.((177|) and (@~9]) yield 

^ = = d 2 i>, (5.3) 

where <9 2 = S^didj. These are the same as the ones 
obtained in GR, and lead to ip — with proper bound- 
ary conditions. It is interesting to note that in GR the 
equation d 2 ip = is obtained from the local Hamiltonian 
constraint, while in the present setup it is obtained from 
the variation of the gauge field A. From this analysis, 
one can see clearly the reason why A is needed in order 
to eliminate the sprin-0 graviton. On the other hand, 
the trace and traceless parts of the dynamical equations 
(|4TT0|) and P~T3l yield 

B = 8 A. (5.4) 

Using the U{1) gauge freedom (I2.3[) . without loss of gen- 
erality, we can set 

SA = 0. (5.5) 

Note that this gauge choice is consistent with our quasi- 
longitudinal gauge tf> = E = [g], because under this 
U(l) gauge transformation, E and (p remain the same, 
as one can see from Eq. (|2.3l) . Then, Eq. (l5.4[) yields B — 
B{x), and the gauge- invariant quantities ^ and $ defined 
by Eq. (|4.3l) are zero, 

* = $ = 0. (5.6) 

Therefore, the scalar perturbations of the metric van- 
ish identically in the Minkowski background. Hence, 
the spin-0 graviton is completely eliminated in the HMT 
setup |33j j - 

It should be noted that Sip is undetermined in the 
present case. However, since (p = 0, it is quite reasonable 
to assume that its linear perturbation also vanishes in 
the Minkowski background. 

VI. CONCLUSIONS 

Recently, Horava and Melby- Thompson [33[ proposed 
a new version of the HL theory of gravity, in which the 
spin-0 graviton, appearing in all the previous versions of 
the HL theory, is eliminated by introducing a Newtonian 
pre-potential <p and a local U(l) gauge field A. Due to 
such an elimination, the dynamical coupling constant A, 
which characterizes the deviation of the kinetic part of 
the action from that of the Einstein-Hilbert, is forced to 
take its relativistic value A = 1. As a result, the theory 
in the IR regime exhibit many features that are quite 
similar to those given in GR. 



Motivated by these remarkable features, in this paper 
we have studied the theory in some detail by assuming 
the presence of matter fields. The potential of the ac- 
tion has been taken to be the one constructed by SVW 
[7J, which represents the most general potential, which 
respects the parity and its highest order of the spatial 
derivatives is six. We have first derived the Hamiltonian 
and super-momentum constraints, given, respectively, by 
Eqs. (|2~TTj) and (|2T3|) . and then the field equations (|27l5)) 
and (|2.16[) . respectively, for the Newtonian pre-potential 
tp and the local U(l) gauge field A. The dynamical 
equations are given by Eq. (|2.18p . while the conservation 
laws of energy and momentum are given, respectively, by 
Eqs.(H211) and (|2~2l3j) . 

Applying the above general formulas to cosmology, we 
have obtained the general modified Friedmann equation 
(I3.8|) and the equation (|3.10l) for the acceleration a. It is 
remarkable that these equations give precisely the con- 
servation law of energy, which takes the same form as 
that given in GR and can be also obtained from the con- 
servation law (|2.24[) . despite the fact that J v and J a are 
non- vanishing, and given, respectively, by Eqs. (|3.12p and 
(I3.13[) . When the spatial curvature is different from zero, 
terms acting as dark radiation and stiff-fluid are present, 
and bouncing universe can be easily constructed from 
these terms. 

We have also studied the scalar perturbations of the 
FRW universe with any given spatial curvature, and 
the linearized Hamiltonian, momentum constraints, the 
equations for the Newtonian pre-potential ip and the 
gauge field A, the trace and traceless parts of the dynam- 
ical equations are given, respectively, by Eqs. (|4.6l) . (14.71) . 
P~g)l , (HU), (ETTU)) , and KW) , while the conservation 
laws of energy and momentum are given, respectively, by 
Eqs.gZEU) and (|4~T51 . 

Applying these formulas to the Minkowski background, 
we have shown explicitly that the metric scalar pertur- 
bations vanish identically, that is, the spin-0 graviton 
appearing in all the previous versions of the HL theory 
is eliminated in the current HMT setup. 

Since the Newtonian pre-potential ip and the gauge 
field A have no contributions to the vector and tensor 
perturbations, the corresponding linear perturbations are 
given precisely by the same equations as those recently 
presented in [ll| in the SVW setup. In particular, in the 
Minkowski background vector perturbations also vanish, 
although the tensor perturbations in general do not . 
These two non- vanishing components represent the mass- 
less spin- 2 gravitons, which are exactly the same as those 
found in GR. 
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